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^ ' Single-particle overlap functions and spectroscopic factors are calculated on the basis of Jastrow- 

type one-body density matrices of open-shell nuclei constructed by using a factor cluster expansion. 
The calculations use the relationship between the overlap functions corresponding to bound states 
of the {A — l)-particle system and the one-body density matrix for the ground state of the A-particle 
system. In this work we extend our previous analyses of reactions on closed-shell nuclei by using the 
resulting overlap functions for the description of the cross sections of (p, d) reactions on the open s-d 
' shell nuclei ^''Mg, ^*Si and ^^S and of "^^8(6, e'p) reaction. The relative role of both shell structure 

and short-range correlations incorporated in the correlation approach on the spectroscopic factors 
and the reaction cross sections is pointed out. 
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OO . I. INTRODUCTION 

o 

I ' Nowadays consistent efforts in the correct treatment of the beyond mean-field nucleon-nucleon (NN) correlations 
, when performing nuclear calculations have been made (see, e.g., [y, |^). While collective phenomena have been known 
' since the early times of nuclear physics, effects produced by short-range correlations (SRC) are more difficult to 
be experimentally singled out. The experimental and theoretical works in the past decade have provided us with 
a much clearer picture on the consequences of these correlations. There are clear signatures of the presence of 
correlation effects on some quantities related to the behaviour of the single nucleon in nuclear medium. For instance, 
^ , NN correlations are responsible for the reduction of the spectroscopic strengths of the nuclear hole states H, ^. 



The theoretical understanding of this fact requires the knowledge of the removal spectral function which can be 
represented in a natural way by both overlap functions and single-nucleon spectroscopic factors . They are directly 
related to observable quantities such as (p, d), (e, e'p) and (7,p) reaction cross sections. The comparison between the 
theoretically calculated and measured cross sections could serve as a test for the proper account of correlations within 
^ the theoretical approach considered. 

Recently, a general procedure has been adopted jsj to extract the bound-state overlap functions and the associated 
spectroscopic factors and separation energies on the base of the ground-state one-body density matrix (OBDM). Ini- 
tially, the procedure has been applied 1^ to a model OBDM 0| accounting for SRC within the low-order approximation 
(LOA) to the Jastrow correlation method B. First, the applicability of the theoretically calculated overlap functions 
has been tested in the description of the ^%)(p, d) pickup reaction |l^ and of the *°Ca(p, d) reaction |9[ |ll|. The 
main result from these investigations is that absolute cross sections of the (p, d) reactions are evaluated without any 
additional normalization in contrast to standard Distorted Wave Born Approximation (DWBA) calculations. Then, 
a detailed study of the cross sections of the electron- and photon-induced knockout reactions on ^^O jl^ and '^''Ca 
PI has been performed. The (e, e'p) and (7,p) reactions are more suitable to test various single-particle (SP) overlap 
functions and NN correlations. The theoretical results obtained in Refs. |l2| reproduce with a fair agreement the 
shape of the experimental cross sections, in particular at large values of the missing momentum, where correlation 
effects are more sizable. Of course, the general success of the above procedure depends strongly on the availability of 
realistic one-body density matrices. 

In the various approaches, the OBDM's are usually constructed for closed-shell nuclei. The ^^O nucleus has been 
studied by variational Monte Carlo and Green function |l^ methods. The treatment of heavier nuclei has not 
yet attained the same degree of accuracy as the light ones. The correlated basis function (CBF) theory, based on the 
Jastrow approach, has recently been extended to medium-heavy doubly closed-shell nuclei using Fermi hypernetted 
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chain integral equations In addition, ^^O and '^''Ca nuclei have been examined by the generator coordinate 
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method |1§. 

There is no systematic study of the one-body density matrix (and related quantities) which includes both closed- 
and open-shell nuclei. For that reason, in expressions for the OBDM's pfr, r') and momentum distributions which 
could be used for both closed- and open-shell nuclei have been obtained. In |20| analytical expressions for the charge 
form factors and densities of the s-p and s-d shell nuclei have been also derived. The p(r, r') was constructed using the 
factor cluster expansion of Clark and co-workers |2^, ^ and Jastrow correlation function which incorporates SRC 
for closed-shell nuclei. This approach was extra pola ted to the case of N=Z open-shell nuclei in Ref. Analyzing 



three different expansions it has been shown in |24 that they lead to almost equivalent results for important nuclear 
properties. The expressions for p(r,r') are functionals of the spherical harmonic oscillator (HO) orbitals and depend 
on the HO parameter and the correlation parameter. The values of the parameters which have been used for the 
closed-shell nuclei ''He, "'^O and ''°Ca were determined in by a fit of the theoretical charge form factor, derived 
with the same cluster expansion, to the experimental one. For the open-shell nuclei ^^C, ^''Mg, ^®Si and '^^S new 
values of these parameters have been obtained to give a better fit to the experimental data Q . Moreover, these nuclei 
were treated as Id shell nuclei and as ld-2s shell nuclei. In the latter case the ^-dependence of the high-momentum 
components of the momentum distributions becomes quite small. 

The aim of the present work is to study the bound-state overlap functions for open-shell nuclei on the basis of 
Jastrow-type one-body density matrices for such nuclei. Next, the resulting overlap functions are tested in the 
description of nuclear observables, namely the cross sections of (p, d) reactions on ^''Mg ^*Si ||2^ and '^^S 
and of the ■^^S(e,e'p) reaction |2^. Such an investigation allows us to examine the relationship between the OBDM 
and the associated overlap functions in that region of nuclei and also to estimate the role of SRC incorporated in the 
correlation approach used in the reaction cross section calculations. 

The paper is organized as follows. A short description of the correlated OBDM of the target nucleus is given in 
Sec. II together with the procedure to extract the SP overlap functions from it. The results of the calculations are 
presented and discussed in Sec. HI. The summary of the present work is given in Sec. IV. 

II. THEORETICAL APPROACH 

A. Correlated one-body density matrix 

In order to evaluate the SP bound-state overlap functions, one needs the ground-state OBDM of the target nucleus 
which is defined by the expression: 

p{r,r')^A j ^'*(r,r2,...,rA)^'(r',r2,...,rA)dr2...drA, (1) 

where ^(ri,r2, ■■■jYa) is the normalized A-nucleon ground-state wave function. The integration in Eq. (|l|) is carried 
out over the radius vectors V2 to va and summation over spin variables is implied. In the present work we start from 
a Jastrow-type trial many-particle wave function 

A 

^{v^,...,VA)^T<^> = Y{f{n,)<^{v^,...,VA), (2) 

where T is a, model operator which introduces SRC, is an uncorrelated (Slater determinant) wave function built 
up from single-particle wave functions which correspond to the occupied states, and f{rij) is the state-independent 
correlation function which was chosen to have the form 

/(r,,) = l-exp[-/3(r,-r,)2]. (3) 

The correlation function goes to 1 for large values of r,j — jr^ — | and to for r^j 0. Apparently, the effect of SRC 
introduced by the function f{rij) becomes larger when the correlation parameter /3 becomes smaller and vice versa. 
For our aims we represent the one-body density matrix p(r, r') by the form 

p(r,r') = ^"^'^^l^l^^'^ = iV(vl/|Orr'|*') = iV(Orr'), (4) 

where ^' — "^{v'^, r^, r^), N is the normalization factor and the integration is carried out over the vectors ri to ta 
and r']^ to r^. The one-body "density operator" Orr' has the form 

A A 

i=l j=jLi 
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The same one-body operator has been used also in ||2^ [Eq. (14) of this reference] for reahstic study of the nuclear 
transparency and the distorted momentum distributions in the semi-inclusive process ^He(e, e'p)X. 

In order to evaluate the correlated one-body density matrix pcor(r,r') we consider firstly the generalized integral 



I{a) ^ (*|exp[a/(0)Orr'] |*'), 
corresponding to the one-body "density operator" Orr' (given by (|5|)), from which wc have 

'din I (a)' 



(Orr') 



da 



(6) 



(7) 



a=0 



For the cluster analysis of Eq. (|^), after considering the sub-product integrals ||2l|, ^ and their factor cluster 
decomposition following the procedure of Ristig, Ter Low, and Clark, one can obtain an expression for the Pcor(r, r') 



p,or{r, r') ~ iV[(Orr')i - O22 (r, r', gi) - 022(r, r', ga) + 022(r, r', gg 



(8) 



The one-body contribution to the OBDM, (Orr')i and the three terms 022(1", r',g;) (/ = 1,2,3), which come from 
the two-body contribution, have the general forms 



(Orr')l =PSD(r,r') = - Vnl{2l + l)(l)*r,i{r)(j)nl{r')Pl{cOSUJrr'), 



(9) 



nl 



and 



022(r,r',a) - 4 ^ TlnihllnjljC^h + l)(2^j + 1) 



h+lj 



4<:;::^;r°(r,r',go - ^ (^,0^,0 1 fco)^<^;;-;f (r,r',a) 



fe=0 



(10) 



where r]ni are the occupation probabilities of the various states, 4'ni(r) are the radial SP wave functions and 
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An 



X / C2;2(^2)0n4/4(^2)exp[-^r2] Ife(2^rr2)r2(ir2 



(11) 



The matrix elements for ^ = 2 and 3 have similar structure and are given in Ref. ||l9|]. In the 

above expression ujrr' is the angle between the vectors r and r' and ik{z) the modified spherical Bessel function. 

Expansion (||) of the present work has one- and two-body terms and the normalization of the wave function is 
preserved by the normalization factor N. The expressions of (Orr')i and the three two-body terms 022(1", r',gi) 
{I = 1, 2, 3) given by Eqs. (g) and (|l^ depend on the SP radial wave functions and so they are suitable to be used 
for analytical calculations with HO orbitals. These expressions were derived for the closed-shell nuclei with N = Z, 
where ry„i is or 1. For the open-shell nuclei (with N = Z) we use the same expressions as in Refs. [|l^ and |^ where 
now < rjni < 1- In this way the systematic study of important quantities such as bound-state overlap functions 
and spectroscomc factors ^ |ll], [T2| can be extended to open-shell nuclei within the factor cluster expansion from 
Refs. [|||l|,@,||. 

It should be noted that a similar expression for Pcor(r, r') was derived by Gaudin et al. |^ in the framework of the 
LOA mentioned in the Introduction. Their expansion contains one- and two-body terms and a part of the three-body 
terms so that the normalization property of the OBDM is fulfilled. 



B. Overlap functions and their relation to the OBDM 



The quantities related to the {A — l)-particle system, such as the overlap functions, the separation energies and the 
spectroscopic factors for its bound states can be fully determined in principle by the one-body density matrix for the 
ground state of the A-particle system ||^ . This unique relationship holds generally for quantum many-body systems 
with sufficiently short-range forces between the particles and is based on the exact representation of the ground-state 
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one-body density matrix g . The latter can be expressed in terms of the overlap functions (pa (r) , which describe the 
residual nucleus as a hole state in the target, in the form: 

p(r,r') = ^C(r)0„(r'). (12) 

a 

In the case of a target nucleus with J'^ = 0+ each of the bound-state overlap functions is characterized by the set 
of quantum numbers a = nlj, with n being the number of the state with a given multipolarity I and a total angular 
momentum j. The asymptotic behaviour of the radial part of the neutron overlap functions for the bound states of 
the {A — l)-system is given by [HI 

(l>7ii]ir) Cnij e^p{-knijr)/r, (13) 

where 

k^i, = ^[2m^iE^,-^'-E^)]'/'. (14) 

In Eq. (|l^) TTin is the neutron mass, Eq is the ground state energy of the target v4-nucleus and E^^^ is the energy 
of the nZj-state of the {A — l)-nucleus. For protons some mathematical complications arise due to an additional long- 
range part of the interaction originating from the Coulomb force, but the general conclusions of the consideration 
remain the same. The asymptotic behaviour of the radial part of the corresponding proton overlap functions reads 

4>nij{r) Cnij cxp[-fc„(jr - r]ln{2knijr)]/r, (15) 



where rj is the Coulomb (or Sommerfeld) parameter and knij contains the mass of the proton. 

The lowest n = no neutron /j-bound-state overlap function is determined by the asymptotic behaviour of the 
associated partial radial contribution of the one-body density matrix py(r, r') {r' — a ^ oo) and Eqs. ( p!2| ) and ( p^ ) 
lead to the expression 

Mr) = r TT' (16) 

6„oij exp(-fc„oij a) /a 

where the constants Cnoij and kn^ij are completely determined by pij{a, a). In this way the separation energy 

and the spectroscopic factor Sngij ~ {4'noij\4>noij) can be determined as well. 

The applicability of this theoretical scheme has been demonstrated in Rcfs. 1^, ||, |ll|, |l^ by means of realistic 
one-body density matrices of ^^O and '^''Ca constructed within different correlation methods. In particular, the 
calculated overlap functions corresponding to the OBDM's from various approaches to the CBF theory jT^, |3|] have 
shown the substantial deviation of their shapes in respect to the Hartree-Fock (HF) wave functions ||l^ . The inclusion 
of short-range and tensor correlations has caused a depletion of the levels below Fermi level which is reflected in the 
substantial reduction of the spectroscopic factors for the quasihole states. 

Thus having the procedure for calculating such important quantities as the overlap functions and the spectroscopic 
factors one can apply it to the one-body density matrices of some open s-d shell nuclei. In order to reveal better the 
nuclear structure properties in this region of nuclei and to draw a parallel with the closed-shell ones, it is desirable to 
test the resulting overlap functions in calculations of one-nucleon pickup and knockout reactions, which is the aim of 
the present paper. 

In the end of this Section we would like to note that the starting point for the present calculations are OBDM's 
which include SRC of central Jastrow type, and which have been computed by a cluster expansion to leading order 
in the correlation function. Similar approach has been very recently used to study the effects of SRC on the (e, e'p) 
exclusive response functions and cross sections |33]. In both calculations the radial dependence of the correlated 



function [Eq. (|^)] is restricted to Gaussian like one. As well known ||10|, |3l|, |3^, 34 , the inclusion of state-dependent 
correlations, in particular of tensor ones, is important when modelling transfer reactions. In addition, long-range 
correlations also should be taken into account in a consistent way. On the other hand, we use in our calculations 
HO single-particle wave functions, which do not have the correct exponential behaviour in the asymptotic region, to 
construct the OBDM. It has been already pointed out in previous works ^ that an important condition of the 
numerical procedure is the exponential asymptotics of the overlap functions at r' — > cxd [see Eqs. dlS) and (fsl)], which 



5 



is related to the correct asymptotics of pijir, r') at r' — > oo. This imposed the development of a method, though not 
unique, in |^ to overcome the difficulties arising from the HO single-particle wave function asymptotics. The condition 
for the correct asymptotics is fulfilled in Q by using of a SP basis obtained with a Woods-Saxon potential. The 
use of HO wave functions in the present work instead of more realistic Woods-Saxon ones following the method from 
allows us, however, to obtain analytical expressions for the OBDM's. At present, no model correlated one-body 
density matrices treating open s-d shell nuclei are available except that of Ref. and, therefore, this work should be 
considered as an initial attempt to study the relative importance of the effects of both SRC and particular structure 
of these nuclei on overlap functions, spectroscopic factors and reaction cross sections. 



III. RESULTS OF CALCULATIONS AND DISCUSSION 

The one-body density matrices (^) constructed by using of factor cluster expansion have been applied to calculate 
neutron and proton overlap functions including NN correlations and related to possible Ip-, Id- and 2s-quasihole 
states in the ^^Mg, ^^Si and '^^S nuclei. We note that in our work the extracted overlap functions are not separated 
in respect to the spin-orbit partners j = I ±1/2. Although the OBDM's do not allow to obtain different results for 
1^3/2 and c?5/2 quasihole states, it is useful in some calculations of reaction cross sections for transitions to these states 
to test the overlap function corresponding to the Id bound state. 

The values of the spectroscopic factors (SF) and of the separation energies for neutrons and protons deduced from 
the calculations are listed in Table 1. The separation energies (|l^) derived from the procedure arc in acceptable 
agreement with the corresponding empirical ones. As a common feature, a substantial reduction of the spectroscopic 
factors of the states which are below the Fermi level (of the independent-particle picture) is observed for all nuclei 
considered due to the short-range NN correlations and their particular open-shell structure. For instance, the values 
of the SF for the states with ^=0 and 1=2 are much smaller than the values obtained for the 1=1 state. The particular 
features of the nuclear structure in the 2s-ld region (in which the Ip state is not a valence one) are responsible for 
the larger reduction of the spectroscopic factors of the levels, though partly occupied, in these open s-d shell nuclei 
in comparison with the "more occupied" Ip quasihole state. It can be also seen from Table 1 that the trend of the 
calculated SF for proton bound states follows that one corresponding to the neutron overlap functions. 

We would like to note the larger reduction of the spectroscopic factors corresponding to quasihole states of the open- 
shell nuclei considered in comparison with the spectroscopic factors for the states in closed-shell ^^O |6[ [l^ and '^"Ca 
1^ nuclei when SRC are taken into account. To discuss this point, firstly we give in Table 1 the spectroscopic factors 
for the Id and 2s quasihole states of the closed-shell '^''Ca deduced from the one-body density matrix calculations 
within the same Jastrow approach fl^ used in the present work. For comparison, the SRC accounted for within the 
LOA to the Jastrow correlation method |^ produce spectroscopic factors of 0.892 for the Id state and 0.956 for the 
2s in '^^Ca nucleus. It turns out that the method from leads to a suppression of the spectroscopic factors for the 
closed-shell '*°Ca nucleus. The reason for the different SF values could be related to the details of the calculations 
with the Jastrow approach when using the two cluster expansions. Particularly they concern the truncation of the 
expansions mentioned already in Section II, which causes a different account for the SRC (although they are of the 
same Jatrow-type) , as well as the different asymptotic behaviour of the OBDM obtained in both expansions. In our 
opinion, the mentioned features of the method from lead to an additional reduction of the values of the SF for 
open-shell nuclei as well. To support this we would like to note the comparison of both cluster expansions given in 
|p4| on the example of the nuclcon momentum distributions n(fc) of various nuclei. The high-momentum tails of n{k) 
obtained using LOA 0| underestimate the corresponding ones obtained with the use of the cluster expansion of Clark 
and co-workers [|l9| for the nuclei considered. In this way the smaller values for the spectroscopic factors derived in 
the present work within for both closed- and open-shell nuclei reflect the stronger accounting for the SRC and can 
be understood. As a common feature, the spectroscopic factors for the s-d open-shell nuclei are smaller than those 
for the closed-shell nuclei. 

Second, in order to check our theoretically calculated spectroscopic factors it is useful to compare them with some 
single-particle occupation probabilities available for the considered ld-2s shell nuclei. The HF calculations performed 
in Q for the {p,p') scattering analyses treating protons and neutrons equivalently yield values of 0.479 and 0.682 
respectively for the 1^5/2 fractional occupancies of ^''^Mg and ^®Si. The 2si/2 proton occupancy of 0.67(9) derived for 
■^^S was reported in Ref. p6t . As can be seen from Table 1 our calculated spectroscopic factors for these states are 
smaller than the resulting occupancies thus satisfying the general property Snij < ^nij^ ^ each Ij subspace the 

spectroscopic factor Snij is smaller than the largest natural occupation number N^j^ ^ . 

Generally, the values of the spectroscopic factors for open-shell nuclei are influenced by the presence of SRC 
accounted for in our OBDM similarly to the case of the closed-shell nuclei. However, in our approach an additional 
reduction of the spectroscopic factors for open ld-2s shell nuclei takes place because the SF contain themselves a 
structure information for the non-complete shell occupancy in these nuclei. This conclusion is confirmed by the check 
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of the normalization condition for the OBDM [Eq. (12)]. 

The DWBA calculations of the pickup reactions on the considered open-shell nuclei were performed using the 
DWUCK4 code assuming zero-range approximation for the p-n interaction inside the deuteron. The same ap- 
proach has been adopted in previous analyses of {p,d) reactions on ^^O |l^ and "^^Ca j|, The values of the 
optical potential parameters have been taken in each case to be the same as in the corresponding standard DWBA 
calculations. 

The standard DWUCK4 procedure is performed by calculating the bound-neutron wave function using the sepa- 
ration energy prescription (SEP) and different sets of proton and deuteron optical model parameters. The optical 
potential is defined to be 

d \ , ^ f h V . Ad 
Vopt^-Vf{x^)~i\W^mD-^ " ' 

where 



Vopt^-Vf{x^)~l[W ^AWD-^]f{xD)~\ Vs.o.{i^.<7)-—J{Xs.o.)+V,, (18) 

V dxD / \ niTiC J r dr 



fix,) = [1 + exp{x,)r\ x,^{r- uA^I^)la,, (19) 

and Vc is the Coulomb potential of a uniformly charged sphere of radius t^A^^'^ ■ The proton and deuteron optical 
model parameters we use are those of Kallne and Fagerstrom ||2^ for ^^Mg and ■^^S and those of Sundberg and Kallne 
||6| for 28Si at £^^=185 MeV incident energy. 

The results for the differential cross sections for the ^''Mg(p, d), ^*Si(p, d) and ■^^S(p, d) reactions at incident proton 
energy Ep—\%h MeV compared with the experimental data are presented in Figs. |l], |^ and ^, respectively. In each 
panel, the results obtained with overlap function, with bound-state wave function following the standard DWBA 
procedure within the SEP and with uncorrelated shell-model wave function are plotted. As can be seen, in general, 
the use of overlap functions derived from the one-body density matrix calculations leads to a qualitative agreement 
with the experimental data reproducing the amplitude of the first maximum and qualitatively the shape of the 
differential cross section. In this case no extra spectroscopic factor is needed, since our overlap functions already 
include the associated spectroscopic factors. In the examples of some reactions the standard DWBA form factor is 
also able to reproduce the shape of cross sections, while the uncorrelated wave functions fail to describe correctly 
neither the size nor the shape of the angular distributions. The differences between the three types of calculations are 
observed mostly at larger angles where the curves corresponding to the uncorrelated case overestimate substantially 
the experimental data. Thus, the relevant importance of SRC included in our approach becomes clear for a correct 
description of pickup processes and this makes it possible to conclude that the replacement of the overlap integral by 
a single-particle wave function is a rough approximation. 

We would like to note that our results for the (p, d) cross sections with the use of overlap functions are obtained 
without any additional normalization while the standard DWBA curves and those corresponding to the uncorrelated 
case need an adjustment by a fitting parameter, i.e., the phenomenological spectroscopic factor. The values of these 
spectroscopic factors are given in Table 2 and can be compared with our theoretically calculated spectroscopic factors 
from Table 1. For instance, the values of the phenomenological spectroscopic factors deduced from the traditional 
DWBA calculations for the transitions to 5/2+ and 3/2+ in ^■^Mg are very close to the empirical values obtained in 
p5| , although a clear physical interpretation could not be done. 

Fig. 1^ shows the differential cross sections for the transitions to the ground 3/2+ state and to the excited 5/2+ 
[at excitation energy E^—'iiAh MeV] and l/2~ [at £'a;=2.76 MeV] states in ^^Mg nucleus. A comparison with the 
experimental data from ]2^ is also made. As can be seen our calculations using the overlap function for the transition to 
the ground 3/2+ state agree fairly well with the experimental angular distribution for the same transition reproducing 
the amplitude of the first maximum and qualitatively the shape of the differential cross section. We emphasize 
that this result is obtained without any additional normalization since the calculated overlap function for 3/2+ state 
already includes the spectroscopic factor of 0.4 (see Table 1). Moreover, at larger angles our calculations lead to 
better agreement with the experimental data comparing with the standard DWBA analysis. In the latter the overlap 
function is replaced by a SP wave function corresponding to a given mean-field potential. In such calculations the NN 
correlations are included approximately by adjusting the mean-field potential parameter values. Another reason for 
the observed discrepancies in |^ is that a spherical potential has been used for generating the bound neutron state 
in DWBA instead of a more realistic deformed one. We would like to mention here the principal role of the overlap 
function for the good description of the differential cross section for the transition to the ground 3/2+ state in ^'^Mg. 
It is obtained on the basis of a correlated OBDM and allows one to account for the short-range correlations in the 
case of pickup reaction. The results of the calculations carried out additionally for the transitions to the excited states 
are less satisfactory for angles larger than 10 degrees. One of the reasons for this is the unrealistic (HO) asymptotic 
behaviour of the Jastrow-type OBDM ||l^ which has been used to calculate the overlap functions following the method 
described in Section II. Concerning the similarity of the shape of the theoretical curves (though different by one order 
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of magnitude) in the upper and middle panels of Fig. || we note that for the ground state the 1^3/2 overlap function 
is needed while for the first excited state the 1^5/2 overlap function is needed. The usage of a common Id overlap 
function in our work obviously gives a larger deviation from the data for the latter case. On the other hand, a large 
fragmentation of the single-particle strengths is a common feature of the studied 2s-lo? nuclei. For example, pickup 
from the lpi/2 an IP3/2 sub-shells in 2s-ld nuclei has been an intriguing topic from experimental point of view 
indicating several Ip states observed in ^^Mg. 

Fig. shows the ground-state differential cross section of the reaction ^*Si(p, c?)^^Si with 1^5/2 neutron pickup, 
while the angular distributions of the reaction ^^S(p, d)^^S for the transition to the ground 1/2+ state and to the 
excited 3/2+ [at E^^IM MeV], 5/2+ [at £^=2.24 MeV] and 1/2" [at E^=7.71 MeV] states in ^^S nucleus are shown 
in Fig. 1^. In general, a quantitative agreement of the calculations with the experimental cross sections in the region 
of the first maximum is obtained for all states of the residual ^^Si and ■^^S nuclei. The presence of admixtures of 
various states observed in the ground state of "^^S [|2^ is responsible for a poorer description of these {p, d) data. The 
behaviour of the theoretically calculated angular distribution and the comparison with the data for the transition to 
the strongly excited l/2~ are similar to those for / = 1 transfer in ^^Mg(p, d)^'^Mg. 

In order to see the sensitivity of the results to the optical potential parameters, we examine the transition to the 
ground 5/2+ state in ^^Si nucleus. In Fig. ^ three different theoretical curves are given in respect to deuteron optical 
potential parameter values used in the calculations. Particularly, the effects of changing the radius of the real part of 
this potential Rd is shown. The best agreement with the experimental data is achieved with the value of Rd— 0.80 fm 
giving also the best DWBA fit in [^. It can be noted that the choice of the radius of the real part of the deuteron 
optical potential within the interval from 0.65 fm to 0.95 fm (considered in the standard DWBA analyses) does not 
influence strongly the good overall agreement of the cross sections obtained by means of the theoretically calculated 
overlap function with the experimental data. Apart from the shown sensitivity of the calculations to the deuteron 
optical potential, in general, we should mention that the {p, d) reaction is more sensitive to the reaction mechanism 
adopted than to the choice of the bound-state wave function. Nevertheless, it is seen from Fig. |4|that our theoretically 
calculated overlap function corresponding to the Id bound state is able to reproduce well the absolute cross section. 

We have analyzed the {p, d) pickup reaction in ^^Mg, ^^Si and ^^S at proton incident energy of 185 MeV. Un- 
fortunately, little is known about the deuteron optical potentials at high energies and the discrepancies with the 
experimental data might be partly attributed to the uncertainties in the determination of the optical potentials. This 
makes also difficult to separate clearly the role of the overlap functions associated with the single-neutron bound 
states. We emphasize, however, the important general necessity the overlap functions to be calculated on the basis of 
realistic OBDM's and only then the sensitivity to all other ingredients of the theoretical schemes to be analyzed. 

It turned out from the previous analyzes of one-nucleon removal reactions pT| , that the quasifree nucleon 
knockout is more suitable to investigate the role of overlap functions as bound-state wave functions. An example 
of electron-induced proton knockout from '^^S for the transition to the ground 2si/2 state of '^^P is given in Fig. 



Calculations have been done with the code DWEEPY |38 , which is based on the nonrelativistic distorted wave 



impulse approximation (DWIA) description of the nucleon knockout process and includes final-state interactions and 
Coulomb distortion of the electron waves ]39[ |. The latter has been treated with a high-energy expansion in inverse 
powers of the electron energy [^ . In the figure the result obtained with the proton overlap function for the 2s state 
of ^^S and the optical potential from [4^ is compared with the NIKHEF data from J^. A reasonable agreement with 
the experimental data for the reduced cross section is obtained. In the analysis of [27| [ the calculations are performed 
within the same DWIA framework and with the same optical potential, but a phenomenological single-particle wave 
fmiction is used with a radius adjusted to the data. We emphasize that in the present work the overlap function is 
theoretically calculated on the basis of the Jastrow-type OBDM of'^^S and does not contain any free parameters. It can 
be seen from Fig. |^ that our spectroscopic factor of 0.5648 gives a good agreement with the size of the experimental 
cross section and, in addition, it is in accordance with the integrated strength for the valence 2si/2 shell in ^^S 
which amounts to 65(7)% of the SP strength obtained using the shell-model bound-state function. The result for the 
•^^8(6, e'p) cross section obtained with the harmonic-oscillator bound-state wave function is also illustrated in Fig. 
It has been computed with the same oscillator parameter value b=2 fm for the 2si^2 ground-state wave function as 
in the original calculations of the OBDM without SRC [|l^ . In order to perform a consistent comparison with the 
result when considering theoretically calculated overlap function, we have applied the same spectroscopic factor of 
0.5648. In this case, the size of the reduced cross section is also reproduced, but the HO wave function gives much 
worse description of the experimental data. Thus, the comparison made in Fig. |^ shows the important role of the 
SRC accounted for in our approach for the correct description of knockout reactions. 

We have performed also calculations for the transition to the excited 1^5/2 state at Ex=2.2?>A MeV and have found a 
good agreement with the shape of the experimental momentum distribution. In our opinion, however, the comparison 
for the Id states is not very meaningful since our approach cannot discriminate between the (^5/2 and (^3/2 states and 
the proton overlap function is calculated for the Id state. From the experimental point of view, it is noted in [ |27t that 
it is also impossible to distinguish strengths originating from both Id shells for the integrated I — 2 strength. In this 
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sense, our spectroscopic factor of 0.6636 fits well the integrated spectroscopic strength in the Id shell up to 24 MeV 
excitation energy (the value 6.04(60) from Table 5 in |^) thus representing a depletion of about 60% of the full Id 
shell. 



IV. CONCLUSIONS 

The results of the present work can be summarized as follows: 

i) Single-particle overlap functions, spectroscopic factors and separation energies are calculated from the Jastrow-type 
one-body density matrices, which were derived using factor cluster expansion, for the ground state of the open-shell 
24Mg, 28Si and ^^S nuclei. 

ii) Taking into account both short-range correlations and specific nuclear structure it is found that the deduced 
spectroscopic factors Snij of the hole states of open-shell nuclei in our particular Jastrow approach are substantially 
smaller than those of the closed-shell ones satisfying at the same time the general relation with the natural occupation 
probabilities Nnij, namely Snij < ^^ij^- 

iii) The absolute values of the differential cross sections of {p, d) reactions on ^''Mg, ^*Si and '^■^S as well as of '^^S(e, e'p) 
reduced cross section are calculated by using the theoretically obtained overlap functions which already contain NN 
correlations. The acceptable agreement with the experimental data shows that this method is applicable also to the 
case of open-shell nuclei. The description of the cross sections might be improved including the deformation effects 
in the even-even Z=N nuclei in the 2s-ld region. 

iv) Using more realistic OBDM's with a correct asymptotic behaviour one can expect that the resulting overlap 
functions will be able to describe more accurately the experimental cross sections of the one-nucleon removal reactions 
on open s-d shell nuclei. In this case more definite conclusions about the role of the SRC in open- and closed-shell 
nuclei can be drawn. 
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Table 1: Neutron and proton spectroscopic factors (SF) and separation energies (e„ and in MeV) calculated on the 
basis of the one-body density matrices |1^ for ^^Mg, ^®Si, '^^S and "^^Ca. Comparison is made with the experimental 
data and Cp^P for the separation energies |6|, p7|, 







neutrons 


protons 


Nucleus 


nl 


en SF 


SF 


24Mg 


Ip 


18.85 19.29 0.9474 


14.13 14.33 0.9798 




Id 


16.94 16.53 0.4026 


11.02 11.69 0.4586 




2s 


18.37 18.89 0.4706 


13.63 14.08 0.5414 


28Si 


Ip 


21.34 21.35 0.8528 


15.21 15.64 0.9788 




Id 


17.29 17.18 0.5071 


11.49 11.59 0.6265 




2s 


17.32 17.96 0.3958 


12.09 12.43 0.4675 


32s 


Ip 


22.54 22.80 0.7454 


15.68 16.29 0.8856 




Id 


16.98 17.33 0.5682 


9.82 10.13 0.6636 




2s 


14.75 15.09 0.4712 


8.53 8.86 0.5648 


40Ca 


Id 


15.45 15.64 0.7691 


8.34 8.33 0.6729 




2s 


17.86 18.19 0.8001 


10.33 10.94 0.8051 



Table 2: Phenomenological spectroscopic factors Sdwba deduced from the standard DWBA calculations and Ssm 
deduced from the calculations with uncorrelated shell-model (SM) OBDM's of Ref. for the reactions considered 
in Figs. |l[|. 



Reaction 


nlj 


Sdwba Ssm 


24Mg(p, d)23Mg 


lc^3/2 


0.25 


0.10 




1^5/2 


2.50 


1.00 




lPl/2 


1.11 


0.75 


28Si(p, d)27Si 


lc?5/2 


2.13 


2.13 


32S(p, d)31S 


2si/2 


0.33 


0.56 




lc^3/2 


1.43 


1.00 




lc^5/2 


4.00 


1.00 




lPl/2 


1.00 


0.11 
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FIG. 1; Diflterential cross-section for the ^''Mg(p, d) reaction at incident proton energy i5p = 185 MeV to the 3/2+ ground 
and to the 5/2+ and l/2~ excited states in ^^Mg. The neutron overlap functions are derived from the OBDM (solid line). 
The calculated standard DWBA curve (dotted line) and the uncorrelated SM result (dashed line) are also presented. The 
experimental data psl are given by the full circles. 



FIG. 2: Differential cross-section for the ^®Si(p, d) reaction at incident proton energy _Ep=185 MeV to the 5/2+ ground state 
in ^^Si. The neutron overlap functions are derived from the OBDM (solid line). The calculated standard DWBA curve (dotted 
line) and the uncorrelated SM result (dashed line) are also presented. The experimental data pq| are given by the full circles. 



FIG. 4: Differential cross-section for the ^*Si(p, d) reaction at incident proton energy i5p=185 MeV to the 5/2+ ground state in 
^'^Si. Line convention referring to calculations using different optical potential parameter values and neutron overlap function 
derived from the OBDM is given (see also the text). The experimental data |2a] are given by the full circles. 



FIG. 5: Reduced cross section of the '^^S(e,e'p) reaction as a function of the missing momentum pm for the transition to the 
1/2+ ground state of ^^P. The proton overlap function is derived from the OBDM (solid line). The result with the uncorrelated 
(HO) wave function is given by dotted line. The experimental data (full circles) are taken from Ref. |27[ . 



FIG. 3: Differential cross-section for the ^^S(p, d) reaction at incident proton energy i5p=185 MeV to the 1/2+ ground and 
to the 3/2+, 5/2+ and l/2~ excited states in "^^S. The neutron overlap functions are derived from the OBDM (solid line). 
The calculated standard DWBA curve (dotted line) and the uncorrelated SM result (dashed line) are also presented. The 
experimental data (Eil are given by the full circles. 
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